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Abstract 

Charm production in pp collisions is considered in the framework of perturbative 
QCD. The values of two parameters, the charm quark mass m c and the QCD 
scale /i 2 , are determined from the comparison of the theoretical calculations with 
experimental data. The RHIC data on charm and beauty production are compared 
with the fcx-factorization approach predictions and with standard NLO QCD. The 
calculated results underestimate the STAR Collaboration data. The role of possible 
nuclear effects is discussed. The predictions for LHC energies are also given. 



PACS numbers: 25.75.Dw, 13.87.Ce, 24.85. +p 



1 Introduction 



The investigation of heavy quark production in high energy collisions is an important 
method for studing the quark-gluon structure of hadrons and the possible nuclear effects 
at early stages of secondary production. The description of hard interactions in hadron 
collisions within the framework of QCD is possible only with the help of some phenomeno- 
logical assumptions which reduce the hadron-hadron interaction to the parton-parton 
one via the formalism of the hadron structure functions. The cross sections of hard 
processes in hadron-hadron interactions can be written in a factorized form [1] as the 
convolutions of the squared matrix elements of the subprocess calculated within the 
framework of QCD with the parton distributions in the colliding hadrons. 

The most popular phenomenological approach is the NLO QCD collinear approxima- 
tion [2-5], where the cross sections of QCD subprocesses are calculated in the Next-to- 
Leading Order (NLO) of a s series. The Fixed Order plus Next-to-Leading-Log (FONLL) [6] 
also resumes large perturbative terms proportional to a™ • \og k (pr /m) with k — n, n — 1, 
where pr and m are the heavy quark transverse momentum and the mass, respectively. 
In these calculations all particles involved are assumed to be on mass shell, carrying only 
longitudinal momenta, and the cross section is averaged over two transverse polarizations 
of the incident gluons. The virtualities q 2 of the initial partons are taken into account 
only through their structure functions. 

The formalism which incorporates the incident parton transverse momenta is refered 
to as the /^-factorization approach [7-10], or the theory of semihard interactions [11- 
19]. Here the Feynman diagrams are calculated taking into account the virtualities and 
all possible polarizations of the incident partons. In the small x domain there is not 
ground to neglect the transverse momenta of gluons, qyr and q2T, when compared to the 
quark mass and to the quark transverse momenta. In such an approach the QCD matrix 
elements of the partonic subprocesses are rather complicated. We have calculated them 
in LO. On the other hand, the multiple emission of soft gluons has been included in our 
calculation. In the /cy-factorization approach the unintegrated gluon distributions are 
used instead of the usual structure functions. 

For the case of heavy ion collisions the linear A-dependence of the heavy flavour pro- 
duction cross sections is usually assumed, so the obtained results are ussually discussed 
in terms of "cross section scaled to pp collision" . This general assumption is experimen- 
tally supported by the result of [20], a ~ A a with a = 1.02±0.03±0.03 (recent results of 
SELEX Collaboration [21] give significantly smaller values of a). In perturbative QCD 
with the factorization approximation [1] the case of heavy ion collisions differs from pp 
collisions only by the change of the usual nucleon structure functions or unintegrated 
gluon distributions by the same functions for bound nucleons. At RHIC energies the dif- 
ference among parton structure functions is known from EMC and NMC experimental 
data [22] and their GLDAP evolution equation analysis [23]. As a result, the nuclear 
effects in QCD for the total cross section of charm production at RHIC are estimated to 
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be on the level of 5-10 % [24] of the linear A-dependence values. 

The experimental data of PHENIX [25, 26] Collaboration obtained at RHIC both 
from pp and nuclear collisions are in reasonable agreement with the NLO QCD and the 
/cy-factorization approach predictions. At the same time, the STAR [27] Collaboration 
data obtained at RHIC in nuclear collisions are in contradiction with standard QCD 
calculations. In the present paper we discuss the accuracy of the theoretical calculations, 
and the plausible reasons for the charm production enhancement on nuclear targets. 

2 Determination of the QCD parameters 

The standard QCD expression for the heavy quark production cross section in a hadron 1- 
hadron 2 collision has the form 

a i2^QQ = f 1 dXa f 1 dxbGa/1 ( Xa ^l) . G b/2 (x b ,n 2 F ) ■ a a ^(s,m Q ,^ R ) , (1) 

Jx a o Jx b0 

where fip is the QCD factorization scale, x a Q = 4toq/s, and x b Q = 4mg/(sa; a ). Here 
G a /i(x a , /Jp) and G b /2(x b , A*f) are the structure functions of partons a and b inside 
hadrons 1 and 2, respectively, and 

a ab ^(s, m Q , /4) = a 2 (/4) • + a^ 2 R ) ■ a™ (2) 

is the cross section of the subprocess ab — > QQ, given by standard QCD as the sum of 
LO and NLO contributions [2, 3, 5]. These contributions depend on the parton center- 
of-mass energy s = {p a + Pb) 2 = x a • x b ■ s, the mass of the produced heavy quark mq 
(actually only on p = Aijiq/s), and the QCD renormalization scale fi 2 R . 

Thus in the QCD calculation of total cross section of charm production three param- 
eters are involved: the heavy quark mass, ttiq, and the two QCD scales, fi 2 F and n R . The 
heavy quark mass rriQ can be estimated as a half of the mass of light quarkonium. More 
accurately, in QCD the masses of heavy quarks should depend on the distance at which 
they are measured. When including this dependence one obtains [28, 29]: 

m c = 1.4 GeV , m b = 4.6 GeV . (3) 

The two QCD scales are usually assumed to be equal in numerical calculations, /i 2 = 
H R = /j, 2 R , and from a general point of view they should be of the order of the maximal 
hardness of the interaction. If the produced heavy quarks are heavy enough, one can 
assume that 

/i 2 ~ itl 2 q , (4) 

though the estimation /i 2 ~ rriq + p\ used in the /c T -factorization approach seems to 
be also reasonable. 
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The next question concerns the uncertainties in the values of these parameters. Usu- 
ally, the accuracy of our knowledge on heavy quark mass is assumed to be about 10% 
[30] (for c-quark), and smaller for 6-quark. The general thinking is that the uncertainty 
of the value of QCD scale /i 2 ~ m 2 Q is about a factor 2 larger. Of course, there is no 
reason why factor 2 should be more reasonable than, say, a factor 1.5 or a factor 3. For 
our purpose, the uncertainties in the values of the discussed parameters should be in any 
case compatible with the values of experimental data including their errorbars. As an 
example, the low boundary for the charm production cross section presented in [31] is in 
total disagreement with all available experimental data, meaning that the uncertainties 
assumed for the values of the parameters were too large. 

We obtain the values of the QCD parameters m c and /i 2 by fitting the results of 
the NLO QCD expression in Eqs. (1) and (2) to the fixed target experimental data in 
[32] and [33], the RHIC point for pp collisions [34], and the UA2 Collaboration point 
[35]. We have used the GRV95 parton distributions [36], which are compatible with the 
most modern analysis (see discussion in [37]). Moreover, the difference to calculations 
performed with other sets of parton structure functions is small in the considered energy 
region. The result of our fit is shown in Fig. 1 by solid curve. 




Figure 1: Total cross section of charm production in pp and pp collisions at ^/s > 25 GeV. 
Experimental fixed target data are taken from [32] and [33], RHIC point from [34], and 
UA2 Collaboration data from [35]. The solid curve shows the result of the fit of QCD 
parameters uiq and /z 2 , and dashed curves show the edges of the uncertainties (the region 
between the y 2 minimal value and the x 2 minimal value plus one). 
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This fit corresponds to the following parameter values: 

m c = 1.48 ± 0.03 GeV , /i 2 = 2.10 ± 0.21 GeV 2 , (5) 

that seem to be very natural. The value of x 2 is 7.3/6 ndf. The uncertainties are shown 
by dashed curves and they correspond to increase the value of x 2 by one from its minimal 
value. 

The same analysis for the total cross section of charm production in up collisions is 
shown in Fig. 2. Here the value of x 2 is 12./ 7 ndf and the parameter values are 

m c = 1.83 ±0.30 GeV , /i 2 = 2.90 ± 0.06 GeV 2 , (6) 

One has to note that in both analysis we obtain /i 2 ~ m 2 , as it is usually assumed 
a priori. Also the two values of m c almost coincide inside the error bars. On the other 
hand, the difference seeming to appear in the /i 2 values can be connected to the not very 
good knowledge of the pion structure functions. 




Figure 2: Total cross section of charm production in -K^p collisions at y/s > 18 GeV. 
Experimental fixed target data are taken from [32] and [33] (they are all shown by 
points). The solid curve shows the result of the fit of QCD uiq and /i 2 , and dashed 
curves show the edges of the uncertainties (the region between the x 2 minimal value and 
the x 2 minimal value plus one). 

3 Charm production in different QCD approaches 

The experimental data on the total cross section of charm production at high energies 
shown in Fig. 1 are presented again in Fig. 3, together with data on interactions with 
nuclei in cosmic rays [38] and RHIC nuclear data [25, 26, 27] scaled to pp interactions. 
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The last cross section data obtained by the STAR Collaboration, as well as cosmic rays 
date, are larger than the theoretical NLO QCD predictions shown by the dash-dotted 
curve. One can also see the different data obtained by the PHENIX and the STAR 
collaborations. This discrepancy is discussed in detail in references [38, 39]. 




Figure 3: Total cross section of charm production in pp and pp collisions at y/s > 25 GeV. 
Experimental fixed target data are taken from [32] and [33], the RHIC point [34], and 
the UA2 Collaboration data from [35] (all are shown by points). PHENIX data [25, 26] 
and STAR data [27] in nuclear collisions are shown by squares and triangles, respectively, 
while cosmic rays data [38] are shown by diamonds. The solid curve shows the theoretical 
result obtained in the /c^-factorization approach, the dashed curve corresponds to NLO 
QCD with only gluon-gluon fusion contribution, and the dash-dotted curve corresponds 
to the total NLO QCD result. 

All other experimental points are in reasonable agreement with the NLO QCD calcula- 
tions performed by taking the values of the parameters presented in Eq. (5) (dash-dotted 
curve), and where the GRV95 parton distributions, compatible with the most modern 
analysis (see discussion in [37]), were used. 

At high energies, the main contribution to the heavy flavour production cross section 
comes from the small x region. Here there is not ground to neglect the transverse mo- 
menta of gluons, qiT and qrr, with respect to the quark mass and to the quark transverse 
momenta, p^- In particular, at high energies and high p iT the main contribution to the 
cross sections comes from the region of qyr ~ Pit, or of q 2 r ~ Pit (see [17] for details). 

For our calculations in the framework of the /cy-factorization approach we have con- 
sidred the heavy quark Q = c,b mass [19] as m c = 1.4 GeV and m b = 4.6 GeV, and we 
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have used QCD scales fi\ = fi 2 F = m\, m T = m 2 Q + p\. All details of the calculations 
in the /^-factorization approach presented below can be found in [19]. The result of 
these calculations are shown by a solid curve in Fig. 3, where one can see that the kr- 
factorization approach gives at not very high energies cross sections which are smaller 
than those obtained by NLO QCD calculations, but which, on the other hand, increase 
slightly faster with energy. 

The cross section of beauty production in pp collisions at RHIC energy y/s = 200 
GeV has been measured by PHENIX Collaboration [43] to be a b i = 3.2l^ IjJ fib, what 
is in agreement with the result of our NLO QCD calculation a h i = 3.06/zb, obtained at 
the same energy y/s = 200 GeV by using GRV95 parton distributions. 

The cross sections of heavy flavour production at RHIC are obtained by the extrap- 
olation of the data measured in the central region to the whole region of rapidities. So 
it seems reasonable to compare the theoretical results to the experimental data before 
extrapolation, what will allow, as a minimum, to avoid a part of the uncertainties. For 
this sake, we present in Tables 1 and 2 the experimental values of da/dy (\y\ = 0.) for 
charm and beauty production, together with the theoretical values obtained both in NLO 
QCD calculations with GRV95 parton distributions and in the fey-factorization approach. 
Again, the STAR data for charm production appear to be several times larger than the 
theoretical results, whereas the PHENIX data for both charm and beauty production 
are in a reasonable agreement with theory. 



Collaboration/ 
Theory 


Reaction 


da/dy (\y\ = 0) 


PHENIX 


Au-Au (m.b.) [40] 
p-p [26] 


143 ± 13 ± 36^b 
123 ± 12 ± 45,ub 


STAR 


d-Au [41] 
Au-Au (m.b.) [41] 
Au-Au (central) [41] 
Cu-Cu (m.b.) [42] 


301 ± 44 ± 67 fib 
267 ± 19 ± 49/xb 
283 ± 12 ± 39/xb 
260 ± 34/ib 


NLO QCD 


p-p 


83^b 


/cT-factorization 


p-p 


56/xb 



Table 1. The experimental and theoretical values of da/dy (\y\ = 0.) for charm produc- 
tion at the RHIC energy y/s = 200 GeV, scaled to pp. 
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Collaboration/ 
Theory 


Reaction 


da/dy (\y\ = 0) 


PHENIX 


P-P [43] 


0.92^;- 


NLO QCD 


P-P 


0.97/ib 


/c-r-factorization 


P-P 


1.62 fib 



Table 2. The experimental and theoretical values of da/dy (\y\ = 0.) for beauty produc- 
tion at the RHIC energy y/s = 200 GeV, scaled to pp. 

In Table 3 we present the predictions for heavy flavour production cross sections and 
midrapidity inclusive densities that we have obtained in the framework of /c T -factorization 
approach. 



Reaction 


Quantity 


yfe = 5.5 TeV 


yfe = 7 TeV 


y/s = 14 TeV 


pp — > cc 


a 


9 mb 


11.3 mb 


21 mb 


da/dy (\y 


= 0) 


1.1 mb 


1.35 mb 


2.3 mb 


pp — > bb 


a 


0.44 mb 


0.57 mb 


1.2 mb 


da/dy (\y 


= 0) 


0.063 mb 


0.079 mb 


1.5 mb 



Table 3. The predictions for charm and beauty production cross sections a, and for the 
inclusive densities da/dy (\y\ = 0.), obtained in the A; T -factorization approach for pp 
collisions at LHC energies. 

4 Transverse momentum distributions 

The STAR Collaboration also presents the p T -distributions of D-mesons produced in 
d+Au collisions and scaled to pp collisions at y/s = 200 GeV. These data, taken from 
[44], are presented in Fig. 4 together with the corresponding theoretical calculations. 

The NLO QCD result shown by the dashed curve is in an evident disagreement (sig- 
nificantly below) with the experimental data on the same level as in Fig. 3. The upper 
curve of FONLL calculations for charmed quark production taken from [44] (dash-dotted 
curve) also underestimates the data. The result of the /c^-factorization approach (solid 
curve) presents a reasonable slope in pr-dependence, but also it underestimates the data. 
In this case, the level of disagreement at high pt is, however, smaller than in the case 
of the total cross sections, what should be connected with the fact that a rather large 
contribution to the total cross section comes from the low-pr region (j>t < m c ), where 
both experimental and theoretical uncertainties are rather large. 
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Figure 4: The STAR Collaboration data for p^-distributions of D-mesons produced 
in d+Au collisions and scaled to pp collisions at ^/s = 200 GeV, together with the 
corresponding calculations in the fey-factorization approach (solid curve), NLO QCD 
(dashed curve), and FONLL [44] (dash-dotted curve). 

It is necessary to note that in Fig. 4 we compare the experimental points for D-meson 
production to theoretical curves for charmed quark distributions. In contradiction to 
what is claimed in [44], on top of fragmentation processes where the momentum of 
the D-meson is smaller than the momentum of the c-quark, also recombination pro- 
cesses where the momentum of D-meson is larger than the momentum of c-quark are 
active. The existence of recombination processes in charm production seems to be ev- 
ident from the experimental data on the asymmetry in yields of the so-called favoured 
and unfavoured D-mesons (see discussions in [45-48]). As a matter of fact, though the 
produced heavy and light quarks have very different transverse momenta, the difference 
in the components of their velocities can be not so large. It is very possible that the frag- 
mentation and recombination processes in charm quark hadronization balance each other 
in the processes with not very high p?, and as a hint of this, the calculated Feynman-rc 
distributions of produced charm quarks in irp collisions are in good agreement with the 
experimental distributions of produced D-mesons (see Fig. 5 in [30]). 

The total cross section of charm production was not measured at Tevatron-collider 
energies. However, data on p T -distributions of D-mesons exist at these energies [49]. 
They are presented in Fig. 5, where one can see that they are in good agreement with 
NLO QCD calculations that use GRV95 parton distributions for charm quarks (dashed 
curve) . 
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Figure 5: The cross sections for .D-meson production in pp collisions at y^s = 1.96 
TeV [49] with \yi\ < 1, together with calculations in the /c^-factorization approach (solid 
curve) and in NLO QCD (dashed curve). 

The solid curve corresponds to the /c T -factorization approach and it slightly over- 
estimates the data, but the agreement should become better in the future when the 
contribution of charmed antibaryons will be added to the experimental data. The re- 
sults of FONLL calculations [44] with fragmentation functions for D-mesons production 
presented in [49], slightly underestimate the yields of D-mesons. Thus, all QCD ap- 
proaches (see also [50]) are in reasonable agreement with the experimental data at yjs 
= 1.96 TeV. Some disagreement of these data with [51] can be expected by the QCD 
parameter values used there, where the charm production cross section increases with 
energy very fast. This should probably result in a rather good agreement of [51] with the 
STAR data presented in Fig. 4, but also in the overestimation of the Tevatron-collider 
data presented in Fig. 5. 

Finally, the theoretical predictions for da/dpr distributions of charm and beauty 
quarks produced in the rapidity window \yi\ < 1 at three LHC energies, ^/s = 5.5 TeV, 
7 TeV, and 14 Tev, calculated in the /c^-factorization approach are presented in Fig. 6. 

5 Nuclear effects in charm production at RHIC 

First, it can be said that the EMC-effect, i.e. the nuclear deformation of parton distri- 
butions [23], should increase the total cross sections for charm production calculated in 
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Figure 6: The predicted cross sections of charm (solid curves) and beauty (dashed curves) 
production in pp collisions at LHC energies, from bottom to top, ^/s = 5.5 TeV, 7 TeV, 
and 14 Tev, in \y\\ < 1, obtained in the /c^-factorization approach. 

the NLO QCD linear approximation at RHIC energy by 5-10 % [24] with respect to the 
linear A-dependence. 

The data of PHENIX and STAR Collaborations scaled to binary interactions are 
presented together in [39]. There one can see the absence of nuclear effects in the 
PHENIX data from pp to central Au-Au collisions, and in STAR data from d-Au to 
central Au-Au collisions, in both cases at the level of 20% accuracy. Concerning the 
PHENIX data presented in [38], they do not present any visible dependence of the 
charm production cross section on the number of collision participant nucleons, N part , in 
Au-Au collisions at different centralities. 

One factor ~2-3 discrepancy [39] exists between the total cross sections of charm 
production obtained by PHENIX and STAR Collaborations. It seems evident that the 
explanation of this discrepancy could be completely connected to some experimental 
problem. However, meanwhile the experimental collaborations will not clarify this point, 
two scenarios can be scripted: either the PHENIX data are right and the NLO QCD 
calculation reasonably describes all experimental data except for two cosmic ray points 
presented in Fig. 3, and then nuclear effects in total cross sections of charm production 
are small, or, on the contrary, the STAR data are correct and the nuclear effects increase 
abot 4-5 times the total cross sections of charm production. In the case the N part 
dependence of these nuclear effects would saturate very fast, even the cosmic ray data 
could be included in the theoretical description. 
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In this second scenario, so large nuclear effects could be connected to large non- 
perturbative contributions (which could be larger than perturbative contributions) to 
high density states, e.g. string fusion [52], percolation [53], strong colour field [54], or 
colour glass condensate effects [55, 56] in the interactions with nuclei. 

In these non-perturbative approaches, the strong colour field inside the cluster formed 
by the overlapping strings produces, above some scale r\ c [53] given by the critical per- 
colation string density, heavy QQ pairs via the Schwinger mechanism in the same way 
as a single string produces light qq pairs. Thus, in the colour glass condensate approach 
[55, 56] the significant scale is the saturation momentum Q s , which grows with energy 
and nuclear size. When Q s > tot(cc), the classical colour field is strong enough and 
it produces heavy pairs QQ. In summary, the production pattern for heavy quarks in 
these non-perturbative approaches becomes similar to that of the light quarks and so 
an overall enhancement of heavy quark production cross section is expected. However, 
to get numerical estimation of these nuclear effects a set of additional assumptions is 
required. 

6 Conclusion 

We show that the /^-factorization approach predictions are in reasonable agreement with 
the NLO QCD calculations for the total cross section of charm and beauty production 
at RHIC energies. 

We obtain a reasonable description of the PHENIX Collaboration data including the 
points obtained from nuclear collisions, and our results are clearly in contradiction with 
the STAR Collaboration data. 

It seems that the main disagreement of our calculations with the STAR Collaboration 
data comes from the low-pr region (p? < 1 GeV/c), where both experimental and 
theoretical uncertainties are rather large. 

The predicted pr-distribution of the produced charm at high pt in the /c^-factorization 
approach is higher than the NLO QCD and FONLL predictions. The /^-factorization 
approach only slightly underestimates the experimental data of the STAR Collaboration. 

Since the NLO QCD calculation basically agrees with the PHENIX data, one would be 
in principle prompted to think that the PHENIX data are right. Either way, the question 
of the existence of a factor ~2-3 discrepancy between the total cross section of charm 
production obtained by PHENIX and by STAR Collaborations has to be experimentally 
unraveled in order to understand the basic mechanisms working in charm production. 

The predictions for LHC energies are also presented. 
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